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Abstract 

We study the effect of the resummation of logarithms for ti production near 
threshold and inclusive electromagnetic decays of heavy quarkonium. This analy- 
sis is complete at next-to-next-to-leading order and includes the full resummation 
of logarithms at next-to-leading-logarithmic accuracy and some partial contribu- 
tions at next-to-next-to-leading logarithmic accuracy. Compared with fixed-order 
computations at next-to-next-to-leading order the scale dependence and conver- 
gence of the perturbative series is greatly improved for both the position of the 
peak and the normalization of the total cross section. Nevertheless, we identify 
a possible source of large scale dependence in the result. At present we estimate 
the remaining theoretical uncertainty of the normalization of the total cross sec- 
tion to be of the order of 10% and for the position of the peak of the order of 
100 MeV. 
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1 Introduction 



The existence of (heavy) quarks with a large mass compared with Aqcd, like 
the top, the bottom, and maybe the charm, makes particularly interesting the 
study of physical processes where a heavy quark pair is created close to threshold, 
because accurate experimental information over the whole near-threshold region 
may allow for a precise determination of some parameters of the Standard Model. 
For instance, the future International Linear Collider (ILC) offers the opportunity 
to study the top quark with unprecedented accuracy [H 13 E] • To fully exploit the 
potential of an ILC in this respect, it is essential that a dedicated measurement of 
the cross section for the production of a top antitop quark pair close to threshold 
is made. Such a threshold scan allows for an extremely precise measurement of 
the top-quark mass and yields information on the top-quark width and the top- 
Higgs Yukawa coupling. The analogous threshold scan for the bb sector is also 
fundamental for non-relativistic sum rules jH IHl El 13 Ej and may lead to accurate 
determinations of the bottom quark mass [Hj. 

The characteristic feature of heavy-quark pair production close to threshold is 
the smallness of the relative velocity v of the heavy quarks in the centre of mass 
frame. This entails a hierarchy of scales fih ^ fJ's ^ f^us where the hard scale fih 
is of the order of the heavy quark mass m, the soft scale /i^ ~ mv is of the order 
of the typical momentum of the heavy quarks and the ultrasoft scale ~ mv"^ 
is of the order of the typical kinetic energy of the heavy quarks. The presence 
of an additional small parameter can be exploited by systematically expanding 
in the strong coupling as and v. Thus, in this context, a next-to-leading order 
(NLO) calculation takes into account all terms that are suppressed by either 
or V relative to the leading-order (LO) result, whereas a next-to-next-to-leading 
order (NNLO) calculation includes all terms suppressed by two powers of the 
small parameter ~ f . The coefficients of this perturbative series contain large 
logarithms logf. In order to improve the reliability of the calculation, these 
logarithms should be resummed. Counting alogf ~ 1 in a NLO or NNLO 
calculation produces a result of next-to-leading logarithmic (NLL) or next-to- 
next-to-leading logarithmic (NNLL) accuracy respectively. There are no leading 
logarithmic (LL) corrections, thus the LO and LL results are the same. 

The expansion as well as the resummation of the logarithms can be organized 
most efficiently by using an effective theory (for a review see Ref. [10,) approach. 
This is done most conveniently by using the threshold expansion ^T] that allows 
to separate the full result of an integral into contributions due to the various 
modes. Denoting the generic integration momentum hj k = (/c°,k), the modes 
that are relevant are: the hard mode ~ k ~ /x/j, the soft mode ~ k ~ 
/is, the potential mode ~ fius, k ~ /x^ and the ultrasoft mode ~ k ~ 
fius- The standard procedure is to first match QCD to non-relativistic QCD 
(NRQCD) ^2] at the hard scale. This corresponds to integrating out the hard 
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modes using the threshold expansion. The resulting theory is then matched to 
potential NRQCD (pNRQCD) O Ej by integrating out the soft modes and 
potential gluon modes. At this stage the theory consists of a non-relativistic 
quark pair interacting through potentials and ultrasoft gluons. 

Within this framework NNLO calculations have been performed by several 
groups (for a review in the case of the top quark see jB]) and quite a few partial 
results needed for a NNNLO calculation have been obtained [1^1 CEl HZl UHl UHl 
|2nil^- Moreover, the existence of the effective theory also provides the necessary 
framework on which to use renormalization group (RG) techniques. These allow 
to resum the large logarithms, logf, that appear as the ratio of the different 
scales appearing in the physical system: log fih/fJ-s and \og Hs/lJ^us- At present, 
the situation is as follows. The computation of the heavy quarkonium spectrum 
is known with NNLL accuracy 1^ (for the hyperfine splitting at NNNLL 
[211 123 ) • Inclusive electromagnetic decays of heavy quarkonium have been com- 
puted to NLL I2B| and for the spin- zero, spin-one ratio to NNLL j2Z|. In 
the case of top-quark pair production a renormalization-group improved (RGI) 
calculation is available |^ I29j. using a somewhat different approach, referred to 
as vNRQCD 30 (in this theory, soft degrees of freedom are kept dynamical and 
the matching from QCD to vNRQCD is carried out directly). However, so far, 
there is no RGI calculation of heavy-quark pair production within the conven- 
tional pNRQCD approach. It is the main purpose of this work to close this gap 
and to provide the ingredients to perform such computations in the case of top 
and bottom-quark pair production. We will report results on tt production near 
threshold, as well as elaborate on the computation of non-relativistic sum rules 
presented in Ref. jHI- Inclusive electromagnetic decays of heavy quarkonium will 
also be considered. These analyses will be complete at NNLO and include the 
complete resummation of logarithms at NLL accuracy and some partial contri- 
butions at NNLL. 

The importance of higher order logarithms can be illustrated in the top quark 
case. In the on-shell scheme, the NNLO corrections turned out to be much larger 
than anticipated and, moreover, made the theoretical prediction very strongly 
scale dependent. If the cross section is expressed in terms of a threshold mass 
l32t l33] rather than the pole mass, the position of its peak is more stable and can 
be predicted with a small theoretical error. This will allow to determine the top 
threshold mass and ultimately the top MS-mass with a very small error. The 
situation is much less favourable regarding the normalization of the cross section. 
The corrections are huge and the scale dependence at NNLO is larger than at 
NLO, indicating that this quantity is not well under control at NNLO. This 
affects how accurately one may obtain the top-quark width and the top-Higgs 
Yukawa coupling. It has been shown that the inclusion of the potentially large 
logf terms is numerically very important and improves the situation regarding 
the normalization of the cross section considerably [2H1 12n]- One of the aims of 
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this paper is to compare with these results. We also investigate to what extent 
the resummation of the logarithms is really required and to what extent the 
improvement we observe in the RGI results is simply due to the partial inclusion of 
higher-order terms. To do so we will produce "NNNLO" results by re-expanding 
full RGI results and dropping terms that are beyond NNNLO. Of course, these 
results are by no means complete at NNNLO, but they contain those NNNLO 
terms that are enhanced by a logarithm. To obtain an estimate of the importance 
of resummation, we then compare these "NNNLO" curves to the full RGI results. 
As we will see, this partial inclusion of NNNLO terms does reproduce the bulk of 
the RGI result except for rather small values of fig- In particular, in the case of 
the top quark the difference between fully resummed and partial NNNLO results 
is small. 

The paper is structured as follows: In Section 2 we present our formalism and 
describe how to perform a RGI calculations in pNRQCD. Some formulae and 
technical details of this section are relegated to the Appendix. In Section 3 we 
apply these results to ti production near threshold. Section 4 contains the appli- 
cation to bottomonium non-relativistic sum rules and inclusive electromagnetic 
decay widths. In the final section, we show our conclusions. 

2 Effective theory 

Within pNRQCD, ip and the fields representing the non-relativistic quark and 
antiquark, interact through a potential and with dynamical ultrasoft gluons. It 
is well known that the leading Coulomb interaction is not suppressed and has to 
be included in the LO Lagrangian which is given by 

where are the colour matrices and the strong coupling is understood to be 
evaluated at the soft scale, = ^^(/is), unless exphcitly indicated otherwise. 
Subleading effects are incorporated in the Lagrangian as corrections to the po- 
tential, 5V , and as interactions of the heavy quarks with ultrasoft gluons. For 
further details we refer to Ref. ^Uj- If we restrict the accuracy of our analysis to 
NNLL, ultrasoft gluons do not appear as physical final states (thought their ef- 
fect is embedded in the RGI running of the matching coefficients of the potentials 
and currents). It follows that the number of particles is conserved (we only have 
one heavy quark and one heavy antiquark) and the problem effectively becomes 
equivalent to do standard quantum mechanics perturbation theory. If we restrict 
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ourselves to study the Hilbert space spanned by the heavy quark- antiquark sys- 
tem in the singlet colour sector we are lead to solve the following equation for 
the associated Green function 

(H{r,p)-E)G(r,r';E)^5{r-r'), (2) 

where 

H^Hc + 5V (3) 

and 

H, = ^-CA, (4) 
m r 

with Cp = {N^ ~ 1)/ (2Ai'c) (where the colour factor Nc = 3). The explicit expres- 
sion of 6V (the correction to the Coulomb potential) will be given afterwards. 
The full Green function G can be solved iteratively in an expansion in the velocity 
by performing multiple insertions of 5V. The LO solution is the Coulomb Green 
function Or and we write 



where 



and 



G(r, r'; E) = Ge(r, r '; E) + 6G{r, r , (5) 



Ge(r,r';£;) = (r|-^— ^lO (6) 



G(r, r'; E) is related to the correlators that appear in the total cross section 
for the production of a heavy quark pair, a{e^e~ QQ) with the centre of mass 
energy = ^/s ~ 2m. This is the key quantity we are interested in. The cross 
section obtains contributions from 7 and Z exchange. In order to simplify the 
discussion we ignore the Z exchange in what follows. The cross section can then 
be written as 

a'^(s)^^-^elR(s), (8) 

where eg is the electric charge of the heavy quark, ckem the electromagnetic 
coupling at the hard scale and the ratio R = a{e^e~ — > QQ)/a{e^e~ — > n'^jJT) 
is expressed as a correlator of two heavy-quark vector currents j^{x) = Q^^Q{x) 

R^s) = ^Im(-^ j rf^a;e'^-^(0|T{/(a;)j^(0)}|0)) . (9) 

In order to compute this correlator, we first express the current in terms of the 
non-relativistic two-component spinor fields ■^^ and x, 

QrQ = cixVv - ^x^^^ i^^f + ■ ■ ■ , (10) 
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where the matching coefficients ci and di are normahzed to 1 at LO. By using 
the equations of motion, Eq. pO|) can also be written in the following way 



hm. V / 



6m 



(11) 



Given the Lagrangian, we solve the corresponding Schrodinger equation and 
ultimately relate the imaginary part of the spin one Green function at the origin 
to R{s) by the equality 



RiE) 



2ATTe%N, 



Q'"^ (cl-c,d,£-^ ImG,=i(0,0;E) 



(12) 



which is valid with NNLL accuracy (where E = ^/s — 2m). We note that in 
evaluating R{E) we expand the expression in Eq. ()12j) and drop all terms that 
are beyond NNLO/NNLL. In particular, we set Ci = 1 in the second term of the 
parenthesis. 

In order to connect with the notation in Ref . JU] , Eq. (fT^ can also be written 

as 



^(^) = Z^S^ImG',=i(0,0;E) 



EM 



Im/ET('^i) + Im^7C('5i 



E 



m 



(13) 



with 



pNR/3 c ^ 
EM I '-'ly 



Im/] 
Im^C(^^i 



"^^qQ^EM 2 



(14) 



2.1 Potential 



The computation of G{0, 0; E) will lead to divergences. These divergences are 
regularized by performing all calculations in momentum space [H^ and using 
dimensional regularisation mD = d + l = 4: — 2e dimensions. Thus, the LO 
Green function at the origin, Eq. (jH)), is understood as 

Ge(r,r';E) = / ^L-^fii^ Ge(p, p'; ^) , (15) 

where G denotes the Fourier transform of the Green function and the insertions, 
Eq. ((7j), are to be evaluated as indicated in Eqs. (jl^ and (jinj) . Using the threshold 
expansion [TT], the Fourier transform of the leading order Green function can be 
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computed as the sum of all ladder diagrams with the exchange of potential gluons 
between the heavy quarks and can be written as 



a(p,p';E) = (27r)'^5W(p-pO^^ \- (16) 



{E - pVm) (p - p')^ {E - p'Vm) 



finite . 



where we have omitted terms that are finite if Eq. ()16|1 is used in Eq. (|T5|l . The 
Green function at the origin has a ultraviolet divergence that manifests itself as a 
pole 1/e. Using MS subtraction and then taking the limit e ^ we find |33ISS1 



a, Cpm^ ( I 1 , -AmE 1 




a(0, 0; E) = - ' f_ ( - + o log — ^ - - + + ^(1 - A) 1 (17) 



where A = Cp dsl {'^y —E /m). 

Higher-order corrections to the Green function can be computed perturba- 
tively in two steps. First, the pNRQCD Lagrangian has to be determined to 
the accuracy needed for the calculation. Second, higher-order corrections are 
computed using quantum mechanics perturbation theory, Eq. 0. Since some 
of the potentials generate singularities in insertions, we have to start from the 
NNLO potential computed in D dimensions |31]. It is essential to manipulate the 
potential consistently in D dimensions, since it allows us to use the same hard 
matching coefficients defined in the MS-scheme obtained in Refs. j36[ IH7j . We 
then bring the potential in a form that is more suitable to be combined with the 
RGI coefficients as presented in Ref. j22]- In particular we use 

^/n|^G.(p,p,)(^ii:ii%)'G,(p3,p.) (18) 

J fJ^{2TxY V(P2-P3) / 

- 2 7r3/2r(l-2e) 

/ n ^ Gc(Pl, P2) I ^^(P3' P4) 

J i=i {^^) rn\p2- P's^^^^ 

to eliminate the Cp term of the non- analytic potential l/g^"*"^"^, present in Ref. |34j . 
The angular momentum operator is generalized to 

(19) 

to be compatible with D dimensional calculations in momentum space. Note that 
an insertion of this operator does not vanish even for an S'-wave (see Eq. 
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but the corresponding contribution could be absorbed into a redefinition of tlie 
matcfiing coefficient of tlie current. 

At NNLL, tlie higher-order corrections to the D-dimensional potential in mo- 
mentum space, SV, can then be written as 

4 

5V = -c,^{27cY5^''\q)-AnCp^ + A7cCF^ (20) 




m 



2 dm^ ^^'^ ' ^^^^ ' 

+ ;j"" |S5.s;]|s^.si] (*--<i!^ 
-5^^^^(|s;.sin-fes|| 

where the colour factor Ca = N^, q = P — p' and contains the corrections to 
the static potential up to NNLL [SHlEni- We will set C4 = 1 due to reparameter- 
ization invariance. We would like to stress that the Wilson coefficients are not 
dimensionless in D 7^ 4. 

The non-relativistic reduction of the spin operators of the potential depends 
on the operators used to single out the physical state we want to study. In our 
case we are using the vector currents, which project to the spin-one state. We 
also have to be careful to use the same conventions than those used to obtain the 
hard piece of the matching coefficients. This produces 0(e) terms multiplying the 

operator. In practise one can do the following replacement for the operator 
with = 2 

\-(<i-m-i)] . (21) 

For S-wave creation there is no contribution from the potentials proportional to 
the S12 and L • S operators. 

The D- dimensional prescriptions used above are irrelevant for the computa- 
tion of the heavy quarkonium mass with NNLL accuracy, which reflects the fact 
that the leading running of the Wilson coefficients is scheme independent. On 
the other hand to specify the D-dimensional prescription of the potential is im- 
portant once they are introduced into divergent potential loops. This is relevant 
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if we want to obtain heavy quarkonium sum rules, or to compute the tt produc- 
tion near threshold with NNLL accuracy, and use computations obtained in other 
places for the hard matching coefficients. We would like to emphasize however 
that we could have chosen a different prescription. This would have changed 
some intermediate-step results but not the physical results. 

Finally, one should note that the potential above is not equal to the potential 
used in Ref. fH3- Nevertheless, they can be exactly related with each other by 
field redefinitions (in four and D dimensions). In particular this means that the 
hard matching coefficients will be the same in both cases. This is actually what 
we expected, since they simply correspond to the effects of the (integrated out) 
hard modes. 

The RGI coefficients of the various potentials are known to the accuracy 
required for evaluating the Green function at NNLL [HHIESI- Note that the strong 
coupling is included in the matching coefficients. Thus D^^^ ~ ^^(l + {as logf )") 

(2) 

and Dx — as{l + (a^ logf )"). For the Coulomb potential ay^, the exact static 



potential at two loop 08] receives additional three-loop LL terms proportional to 
alC\{as log v)" [Sn]. The explicit form of all matching coefficients can be found in 
Ref. (221 • However, we have changed the basis of potentials compared to Ref. 
This affects the matching coefficient -D^^], which now reads {fius = f^l/f^h) 



Dd,si^^us) = (2 + Cd(^.)) + — 4s(^s) + 7^4s(/i. 

o on \ Up 



(22) 



For the other potentials, the expressions obtained for their Wilson coefficients in 
Ref. |39t l^ hold, where one may also find the expressions for the RGI coefficients 
(cd, d^s, dss, . . .) of the NRQCD operators. We repeat them in the Appendix for 
ease of reference. 

In this paper we also include an additional contribution in 5V , the electro- 
magnetic Coulomb term 

5V -^5V . 23 

This will give rise to a NLO term ~ a-Eu/v from single potential photon exchange 
and a NNLO term from double potential photon exchange. Strictly speaking, the 
coupling in Eq. should be evaluated at the soft scale and not at the hard 
scale, but this effect is beyond NNLL. 

This completes all the corrections needed at the Lagrangian level for the 
computation of the imaginary part of the Green function with NNLL accuracy. 
What is left is to obtain the RGI expressions for Cs and rf^, where s = 0, 1 labels 
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the spin. The matching coefficient Cs is needed at NNLL whereas dg is only 
needed at LL. 



2.2 Direct ultrasoft effects to c. and d. 



Most of the uhrasoft contribution to the running of Cg and dg comes in a indirect 
way, through the running of the potentials in the anomalous dimensions of the 
RG equation. Nevertheless, there are some genuine ultrasoft effects that have not 
been considered so far. They are due to the appearance of some energy dependent 
potentials with the structure 



(24) 



where Z is the normalization correction to the heavy quarkonium propagator. 
Therefore, Z^^'^ corresponds to the normalization of the field that represents the 
heavy quarkonium. SV^s was not included in Eq. (j^UI) . as there only energy 
independent potentials were considered. Its effects could be reabsorbed in field 
redefinitions of the fields that represent the heavy quarkonium. Therefore, they 
have no consequences in the spectrum. Nevertheless, these field redefinitions 
change the vertex interaction of the heavy quarkonium with photons producing 
changes in Cg and dg and they have to be considered in our computation. The 
leading logarithmic corrections to Z^^"^ were obtained in Eqs. (15-18) of Ref. J7] 
(for some partial results see also Ref. (Ifij). They produce corrections of order 
ttglogOs. We can obtain the RGI expressions for them and generate terms of 
order log'-^^"'' ag. Note that the running goes up to the soft scale because 

one gets expressions of the form (r = 0| logr|p) and the logarithm gets the scale 
of p. The corrections to the Green function due to Eq. (j^^ read 



SGi 



6Go = — ^lo: 



SC. 



4 



log 



as(/is) 
as(/i«s) 



ClC 



H,-E 4 
1 1 



mr He — E 



3 ^ 



CfCa 



C, 



The correction coming from (5G3 can be included in dg. 



dg dg 



16C; 



log 



as (/is 



(25) 
(26) 
(27) 

(28) 



The correction from is zero. Finally, the correction from bG\ could be ab 

as(/is! 



sorbed in Cg. 



as (/it 



(29) 
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Note also that these changes are equivalent to a change in the NNLO and LO 
anomalous dimension of the RG equation describing the running of and ds 
respectively. 



2.3 Running of ds 

The LL running of Im^fg^JT can be obtained in two steps. In the first step one 
computes its soft running. This has been done in Ref. ^OlE]- For the explicit 
result see Eqs. (C.18,C.19) in Ref. [Hj. The ultrasoft running can be obtained 
from Eq. ()28|) . Adding everything together the LL running of Im^f^^^ reads 



16 



Imc/gM^(^5o)(/i/,) - — CFlm/|M^(^S'o)(/i/,)log 

6pQ 



Im^gM^(^5i)(/ife) - -^C^Im/^M^(^5i)(/i/,)log 



«s(/ift) 
«s(yU/t) 



(30) 



(31) 



The last result agrees with the LL result obtained in Ref. |29] . For the matching 
coefficient di as defined in Eq. (fTUj). this entails 



diifius) = IH n — log 

Po 



(32) 



2.4 Running of q 



The running of Cg is not yet known with NNLL accuracy. It is dictated by the 
solution of the RG equation (the LO anomalous dimension is zero) 



d 
dfi, 



NLO I ,,NNLO 



logc. = 7r" + 7, 



+ 



The structure of the solution reads 



(33) 



(34) 



Expressions for Cs{fih) at two loops in the MS can be found in Ref. [HE] for s = 1 
and (almost complete) in Ref. [SZI for s = 0. 

The expression for the basis of potentials used in this paper reads 



NLO 



CI 



a, + ( 2 - -s{s + 1) ) Z^g, - ^Dli + AD\ 



)(2) 



CaCf 



(35) 
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The expression of T^^^ is known |2H1 123 • For r^'^^^ only the spin-dependent 
term is completely known |2Z]. Its contribution to the spin-one case reads 



■^NNLL/ \ _ pNNLL , ^1^2 O^sjf^s) ^,(2) / n 11^2 



where T^J^'"^'"'" corresponds to the result quoted in Ref. ^7\. Note that we use a 
different expression for the spin- dependent term than just F^J^^'"'". The expres- 
sion above corresponds to the spin- dependent contribution to the vector current 
matching coefficient in the MS scheme. Numerically this contribution is small 
compared with others. 

Besides the spin-dependent correction, we have also incorporated the follow- 
ing spin-independent corrections at NNLL order (to these corrections one obvi- 
ously has to subtract the spin-dependent piece that has already been included in 

r-pNNLLv 
ci,SD)- 

a) those that appear from the exponentiation of the NLL term and formally are 
NNLL (see Eq. (jBl), 

b) Effects due to the two-loop beta running of Og. They produce the following 
correction: 

.NNLL/ X A r^^'''^ ^«s,NLO 



«s^(/^.)5r--(;.,) = f^, I ^ (37) 

_27r da,^ da,{^il/^iH) ^7^° 

(3q JasM as ^ d(3i das{iJ.l/nh) ' 

The last term is generated from the fact that in the determination of r^^^{fis), 
the relation 

«s(/x.) ^(2-zA)' 

where 2;^° = as(/Us)/as(/^/i), was used. This relation is only true at one loop and 
has to be corrected by jSi terms if a NNLL accuracy is demanded. The numerical 
impact of these corrections is small. 

c) We have also incorporated the corrections proportional to oi, the one- loop log- 
independent term, that appears in ay^ — as(/^s) (l + ^^i)- These corrections 
can be deduced from the computation of the NLO anomalous dimension. They 
read 

^2/„ UrNNLL/ N 1^ P^^'^ C?"s / NLO , C^CV ^(1) ^2 "s \ /oqN 

o^Af^hjdT (/ij = -7-5-ai I ^ — lie, + ^^^s ' - Cp-r ■ (39) 

The numerical impact of these corrections is small. 

Finally, the inclusion of the electromagnetic corrections produces some cor- 
rections to C\. Counting oem ~ the one-loop exchange of a hard photon 



11 



contributes at NNLO and is taken into account by 

2e2 



cM-^cM---^^ (40) 

TT 



for the spin-one case, whereas for s = we have 



eoW-eoW-^ l-^r . (41) 



TT 



2.5 Green function 

Once the RGI potential and current matching coefficients are available, we are in a 
position to use standard quantum mechanics perturbation theory to compute the 
higher-order corrections to Im[G(0, 0; E)] via insertions of the potentials. This 
calculation has been done in momentum space using dimensional regularization. 
For the terms suppressed by two powers of ~ f in 6V, Eq- ipUj) . it is sufficient 
to consider a single insertion, 

6GiO, 0;E) = Jl[ ^ Ge(pi, P2, E) 5V(p2, Ps) ^(Pa, P4; E) , (42) 

whereas for terms suppressed by only a single power of ~ f we have to compute 
double insertions as well, 

5G(0,0;E)= (43) 

/ n TTT^ G'c(pi, P2; E) 5y(p2, Ps) ^^(ps, P4; E) 5y(p4, Ps) ^^(ps, Pe; E) . 
•' i=i l^^J 

Since the RGI does not alter the structure of the interaction terms, the insertions 
can be taken directly from Ref. In Appendix^ we list the integrals that are 
needed and present the explicit results obtained in Ref. [23 . 

In the calculation described so far, the on-shell scheme for the heavy quark 
mass has been implicitly assumed. In order to avoid the bad convergence be- 
haviour inherent to this scheme, it is necessary to rewrite the expressions in 
terms of a threshold mass which is free of the renormalon ambiguity. In this 
paper we will consider the cases of the potential subtracted (PS) mass and 
renormalon subtracted (RS) mass [32] ■ We will use the difference between both 
schemes as an indication of the scheme dependence of our results. 



3 The case of the top quark 

In this section we apply the previous results to the case of top-quark pair produc- 
tion near threshold at a future linear collider. Due to the large width Tt of the 
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top quark there are no bound states and the toponium resonances are smeared, 
resulting in a smooth curve for the cross section with a broad peak as the rem- 
nant of the would-be 15* bound state. Using perturbation theory we can reliably 
compute the cross section as a function of the energy p. This may lead to accu- 
rate determinations of the top mass and its total decay width by measuring the 
position and normalization of the peak. 



The tt pair will be dominantly produced via e~^e 
production cross section may be written as [S] 



7^ 



EM 



3s 



+ F''(s)R''(s] 



tt. The total 



(44) 



where 



R^^s) = ^Im J d'xe"^-'{0\T{j;{x)f>'m\0) 
R^{s) = ^Im (^-i J rf^xe^''-^(0|r{j;(x)/'^(0)}|0)^ , 



(45) 



and {j'^) is the vector (axial-vector) current that produces a quark-antiquark 
pair defined by Eq. (fTnjl . With both 7 and Z exchange the prefactors in Eq. 



are 



where 



F"(sl 



T/ - 2ef 



2 _ 2s VeVqeq s"^ [yl + al)vl 



s^ {vl + al)a 



m 



2 \2 

z) 



sin^ e 



w 



(46) 



(47) 



2 sin 9w cos 9w 2 sin 9w cos 9w 

Here e/ is the charge for fermion /, is the third component of weak isospin, 
9w is the weak mixing angle, and mz the mass of the Z. Here we will focus 
on i?'^, since it gives the dominant contribution and we are mainly interested in 
studying the impact of the resummation of logarithms on the convergence and 
scale dependence of the perturbative expansion of R". A study of i?" can be 
found in Ref. 

Since Tt ~ maEu ~ ''tT'v'^ and the propagator of a potential heavy quark 
scales as f ~^ the effects due to the width of the top quark are LO effects and 
have to be taken into account by modifying the propagator E — p^/ (2m) — 
E + iTt — p^/ (2m). This amounts to replacing E E + iTt in Eq. (^. As noted 
in Ref. [Tj, this is only correct at LO. Higher-order electroweak corrections have a 
much richer structure and it is not possible any longer to formulate the problem 
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in terms of a ti final state. In particular, at NNLO there are interference effects 
(between double and single resonant processes), QED radiation effects and non- 
factorizable corrections (non-trivial interconnections between the decay products 
of the top quarks with the remainder of the process). Any consistent approach 
beyond LO has to introduce additional operators with fields corresponding to the 
incoming electrons and the decay products of the top quarks, and link higher- 
order corrections of the ip^ip and x'^X operators to three-point and higher-point 
vertices. Even though there is an effective theory framework available for sys- 
tematically taking into account these corrections |12], a full explicit calculation 
of electroweak effects at NNLO is still lacking. If we are interested in the total 
cross section only, the situation is somewhat simpler, since the non-factorizable 
corrections cancel j3H] , and the replacement E ^ E + iTf becomes correct with 
NLO accuracy. Furthermore, some of the electroweak corrections have be taken 
into account by including them in the matching coefficients |44j. In this article 
we restrict ourselves to the usual shift E ^ E + iVt- 

In Figure HI the fixed-order results are compared to the RGI results obtained 
using the procedure described above and using the PS mass with the subtraction 
scale fip required for the definition of the threshold mass set to 20 GeV. The 
plots were produced with an energy dependent soft scale = Am\J E'^ + F^. 
In order to obtain a first rough estimate of the theoretical uncertainty we show 
the cross sections as bands obtained by variation of the soft scale in the region 
30 GeV < /is < 80 GeV, where these numbers refer to the scale Hs sX E = 0. 
For the hard and ultrasoft scale we take our default values Hh = fnps = 175 GeV 
and fius = yUs/yU/i- Our results are qualitative consistent, but not equal, to those 
obtained in Refs. [23 121] • This agreement is not trivial, since the ingredients 
included in the NNLL analysis are different in ours and their computations. This 
may indicate that, even if the NNLL evaluation is incomplete, the qualitative 
features will hold in the complete result. We observe that the scale dependence 
is much reduced once the logarithms are taken into account and reduces from 
LL to NLL to NNLL. We also note that the size of the corrections decreases for 
the RGI results, in particular the NNLL band is much closer to the NLL band. 
However, the NNLL band does not overlap with the NLL band, indicating that a 
theoretical error estimate relying on the scale dependence alone is too optimistic. 
Therefore, we consider other possible source of errors in what follows. 

The variation of the soft scale has been stopped at fig = 30 GeV which might 
seem to be a rather large value. In fact, the dependence of the normalization of 
the peak as a function of fig is very smooth at NNLL, up to a value of fig ~ 25 GeV, 
where it abruptly changes^. At NLL, already for fig ~ 30 GeV there is a rather 
large variation of the normalization. This is illustrated in Figure |21 which shows 
the fig dependence of the normalization of the peak at LO/LL, NLO, NLL, NNLO 

^As we will see a similar pattern also appears for the bottomonium decays. 
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Figure 1: Threshold scan for ti using the PS mass, mps(20GeV) = 175 GeV. The 
upper panel shows the fixed order results, LO, NLO and NNLO, whereas in the lower 
panel the RGI results LL, NLL and NNLL are displayed. The soft scale is varied from 
Hs=30 GeV to /is=80 GeV. 



and NNLL. As observed in Ref. [20] the situation for small soft scales may be 
remedied if multiple insertions of the Coulomb potential are taken into account. 
Since in our result these, formally, higher-order contributions are not taken into 
account, we refrain from using scales /i^ < 30 GeV and take /i^ = 40 GeV as 
our default value for the soft scale unless stated otherwise. The "NNNLO" result 
also depicted in Figure El is obtained by re-expanding the full RGI result and 
keeping only terms that are NNNLO, but dropping those of even higher order. 
The difference between this result and the full RGI result is small except for 
scales fig < 30 GeV, which is outside the range we use. 
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Figure 2: The normalization of the peak of the RGI threshold cross section as a 
function of the soft scale fj,s ■ The vertical dashed lines show the limits of variation used 
in Figure n 




Figure 3: The position of the peak of the RGI threshold cross section as a function of 
the soft scale /x^- The vertical dashed lines show the limits of variation used in Figure ^ 



We also present a similar plot to show the scale dependence of the position 
of the peak in Figure 01 From this plot we can see that the NNLL resummation 
of logarithms significantly improves over the fixed-order NNLO evaluation. The 
scale dependence is reduced and also the convergence of the perturbative series 
is better, which is reflected in a smaller difference between the NLL and NNLL 
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result versus the NLO and NNLO result (this also happens for the normalization 
of the total cross section). The difference between the "NNNLO" and NNLL curve 
is again very small for fis > 30 GeV. From this plot we estimate the theoretical 
error for the determination of the position of the peak (which is related to the 
determination of the top mass) to be of the order of 100 MeV. 

In the fixed-order calculations the error is usually estimated from the soft- 
scale dependence, even though one may think that the dominant uncertainty 
came from the magnitude of the correction. This is potentially more of a problem 
after the resummation of logarithms, since the soft scale dependence is much less 
severe. Therefore, we have to be more careful with other sources of uncertainties. 
In particular, the missing ultrasoft contributions make it important to consider 
the dependence on the other scales as well. Since they are correlated, we can 
only vary /i/j together with and, thus, consider in Figure El the dependence 
on the hard scale fih, setting /i^ = 40 GeV. Variation of the hard scale around 
its natural value fj,h = m hj choosing 100 GeV < fih ^ 250 GeV results in a 
scale dependence that is considerably larger than the soft scale dependence. We 
note that this error is compatible with the magnitude of the difference between 
the NLL and NNLL result. It is to be expected that the situation improves once 
all ultrasoft logarithms at NNLL are taken into account, but at this stage the 
rather large dependence of the cross section on fi^ has to be taken into account 
if a theoretical error is assigned. We also note that the dependence on fih only 
enters at NLL, thus the LL "band" in Figure E] is simply a line. 
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Figure 4: Dependence of the ti threshold scan on the hard scale ///i, using the PS mass. 
At NNLL (NLL) the lower (upper) curve corresponds to /i/, = 250 GeV, whereas the 
upper (lower) curve corresponds to [i^ = 100 GeV. 

Finally we turn to Figure where we display the effects due to the QED 
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Figure 5: Effects of the QED corrections to the tt threshold scan. The hard ans soft 
sales are chosen as = mps = 175 GeV and Hs = 40 GeV. 



corrections with the default choice for the scales, fih = and fig = 40 GeV. As 
previously mentioned, the QED effects enter at NLL (thus there is no effect on the 
LL curve) and, compared to the desired accuracy (top quark mass measurement 
with an error 5m < 100 MeV), they are large. They change the normalization 
by up to 10% and result in a shift in the extracted MS mass of up to 100 MeV, 
making it mandatory to include them. We note that we have not changed the 
definition of the PS mass. Thus the shifts shown in Figure El are physical effects 
and are not compensated if the PS mass is related to the MS mass. 

We end this section by mentioning that all results presented here can be 
worked out in different threshold mass schemes. The qualitative features of the 
results obtained are similar. We illustrate this point using the RS and RS' mass, 
again setting fip = 20 GeV, and depict the threshold scan in Figure IHl Those 
plots, together with the plot in the PS scheme (lower panel in Figure Q), allow us 
to visualize the effect of using different threshold masses. In principle, the value 
of the threshold masses can run from being numerically close to the MS mass (but 
not too close otherwise power counting is broken) to being numerically close to the 
pole mass (but again not too close, otherwise the renormalon cancellation does not 
take place). Within this allowed range, the threshold mass definitions numerically 
closer to the MS mass have a smaller scale dependence. The price paid is that 
the magnitude of the corrections (between different order in perturbation theory) 
is larger. For the specific threshold masses we use, the RS mass is the one closest 
to the MS mass whereas the RS' mass is closest to the pole mass. The PS mass 
is in the middle, though somewhat closer to the RS mass. Once experimental 
data is available it will be useful to consider different mass definitions to obtain 
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Figure 6: Threshold scan for tt using the RS and RS' mass, mRs(20GeV) = 175 GeV, 
mRS'(20 GeV) = 175 GeV. The upper panel shows the RGI results LL, NLL and NNLL 
with the RS and the lower panel with the RS' mass. The soft scale is varied from 
Hs=30 GeV to /is=80 GeV. 



an estimate of the corresponding error. 



4 The case of the bottom quark 

We can also apply the results of Section 2 to a variety of observables for bb systems. 
Likely, the theoretically cleanest observable on which one can use our results is 
non-relativistic sum rules, which has already been considered in Ref. 0, where an 
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accurate determination of the bottom mass was obtained. Here we will consider 
inclusive electromagnetic decays (either to e^e~ or to 77) of the bottomonium 
ground state for which we will provide analytic formulae. These can be obtained 
from the results obtained for the non-relativistic Green function. The spin-one 
decay has the following structure 



r(T(n5) ^ e+e") 



iOTT 

3 
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The corrections to the wave function at the origin are obtained by taking the 
residue of the Green function at the position of the poles 

0ri)(O)|' = |0W(O)f (l + 5<^(r=i)) = Res a=i(0,0;i?), (49) 

E — En 

where the LO wave function is given by 
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The corrections to produced by 5V have already been calculated with 

NNLO accuracy |S1 IHj in the direct matching scheme. One can also obtain them 
in the dimensional regularized MS scheme with NNLL accuracy by incorporating 
the RGI matching coefficients. One obtains then the following correction to the 
wave function 
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For completeness we also give the corrections to the wave function induced 
by the QED effect, Eq. They read 
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Of course, the corrections to the matching coefficients of the current, Eqs. (j4(Jj) 
and (j4H) have to be taken into account as well. Due to the additional suppression 
Cq = 1/9 these QED corrections are numerically not very important. 

From these expressions one can obtain the wave-function correction for the 
spin zero case 

50(r=°) = 50^'^ + ^<\>n , (54) 

where 
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7 



n 



(55) 



by using the results from Ref. j^j. Therefore, the decay of the pseudoscalar 
heavy quarkonium to two photons reads (c/q = di) 
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We now perform a phenomenological analysis of these results. We restrict 
our analysis to the ground state of bottomonium. For the mass of the r)b{lS), 
we use the T(IS') mass, which is consistent to the order of interest. When we 
perform the numerical analysis, we expand the expressions (except for the overall 
factor l/A^^^(i5)/^j^(i5)) in Eqs. PH|) and (jKHjl . and drop sub leading corrections, in 
order to work strictly at LL, NIL (NLO) and NNLL (NNLO). In particular, this 
applies to the relativistic correction proportional to d^, where MY(is)/r}tiis) ~ 2"^ 
is replaced by Ei^^ = —mCpag/4. The expressions above have been written 
in the pole scheme. Therefore, we change to a threshold scheme suitable for 
the renormalon cancellation. We consider three possible schemes: the PS, the 
RS and the RS' scheme with the subtraction scale fip = 2 GeV. The numerical 
values of the bottom-quark mass to be used in the various schemes are take from 
Ref. and are given by mps(2 GeV) = 4.515 GeV, mRs(2GeV) = 4.370 GeV, 
and mRS'(2 GeV) = 4.750 GeV. Apart from the different numerical values for m 
to be used, the first change in the formulae appears (at most) at NNLL (NNLO) 
and is due to the shift 
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where mx represents a generic threshold mass, F'-"-' is the decay width at lowest 
order and the shift in the mass 6mx ~ mv'^. We notice that in the RS' scheme 
Smus' = at order mv'^. Therefore, in this scheme, with the current precision 
of the calculation, the expressions for F are equivalent to those in the pole mass 
scheme. 
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Figure 7: Prediction for the T{1S) decay rate to e+e . We work in the RS' scheme. 
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The results for the vector and pseudoscalar decay can be found in Figure [7| 
and Figure IHl respectively. From the numerical analysis, we find that the NNLL 
corrections are huge, especially for the rjb^lS) 77 decay. The result we obtain 
for this decay is compatible with the number obtained in Ref. P^. This is some- 
what reassuring, since in that reference the ratio of the spin-one spin-zero decay 
was considered, which was much more scale independent, as well as more conver- 
gent (yet still large) than for each of the decays themselves. This agreement can 
be traced back to the fact that, for the spin-one decay, for which we can compare 
with experiment, we find that the NNLL result improves the agreement with the 
data. Overall, the resummation of logarithms always significantly improves over 
the NNLO result, the scale dependence greatly improves, as well as the conver- 
gence of the series. On the other hand the problem of lack of convergence of the 
perturbative series is not really solved by the resummation of logarithms and it 
remains as an open issue. Due to the lack of convergence we refrain from giving 
numbers (and assigning errors) for our analysis. In this respect we can not avoid 
to mention that, whereas the perturbative series in non-relativistic sum rules is 
sign-alternating, is not sign-alternating for the electromagnetic decays. Finally, 
we would also like to remark the strong scale dependence that we observe at 
low scales, which we believe to have the same origin than the one observed in ti 
production near threshold in the previous section. 

As for the tt case the "NNNLO" curve, which contains the logarithmically 
enhanced NNNLO terms, follows closely the RGI curve up to a certain value of ns- 
For smaller soft scales, the two curves start to deviate. However, in the bb case, 
this deviation takes place already for reasonably large values of fig ~ 2.5 GeV. 
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This seems to suggest that resummation is rather more important in the bb case. 
However, we have to keep in mind that some logarithms at NNLL are still missing 
in our analysis and this might well affect this conclusion. 

We also illustrate the dependence on the threshold mass evaluation in Fig- 
ure El Compared with the uncertainty due to the lack of convergence of the 
perturbative series, the dependence on the threshold mass is negligible. 

1.75 




Figure 9: Prediction for the T(15) decay rate to e+e~ at LL, NLL and NNLL for the 
PS, RS and RS' mass. 



5 Conclusions 

We have studied the effect of the resummation of logarithms for ti production 
near threshold and inclusive electromagnetic decays of heavy quarkonium. This 
analysis is complete at NNLO and includes the full resummation of logarithms 
at NLL accuracy and some partial contributions at NNLL accuracy. 

Compared with fixed-order computations the scale dependence and conver- 
gence of the perturbative series is greatly improved for both the position of the 
peak and the normalization of the total cross section of ti production near thresh- 
old. Nevertheless, we identify a possible source of a large scale dependence in the 
result. While the result is very stable with respect to the variation of the soft 
scale fig it shows a considerably larger dependence on the hard scale /i^. This 
might well be related to the fact that the variation with respect to fih is corre- 
lated with the ultrasoft scale variation and the ultrasoft logarithms are not fully 
included at the NNLL level in our result. At present we estimate the remaining 
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theoretical uncertainty of the normahzation of the total cross section to be of the 
order of 10% and for the position of the peak of the order of 100 MeV, based 
on the difference between the NLL and NNLL plots, as well as on the hard scale 
dependence. We note that this estimate of the theoretical error is somewhat 
larger than in Ref. jlSj, due to the fact that we use more possible sources of the 
theoretical error. We would like to remark that the use of the RGI provides us 
also with a significantly improved determination of the top mass. 

For the inclusive electromagnetic bottomonium decays the corrections are 
even larger than in the top case. In the case of the T(15') —>■ e~^e~ , they bring 
the final prediction into better agreement with experiment. For the rjiy — > 77, 
there is no experimental data at present and the presence of huge corrections 
make a reliable theoretical prediction difficult. 

Finally we remark that the NNNLO part of the RGI results is numerically 
considerably more important than the even higher-order terms, in particular for 
the top case. Thus it would be highly desirable to combine a full NNNLO eval- 
uation with a complete NNLL result in order to obtain a satisfactory theoretical 
prediction of heavy-quark pair production near threshold. 
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A Insertions 



In this appendix we discuss how to obtain the perturbative expansion of the 
Green function. For the loop integration measure we use 



dpi = 
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such that the MS-scheme corresponds to minimally subtracting the poles in e. 
We start by writing the potential, Eq. (|^. for the case of a spin triplet S wave 
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where we have used the exphcit expression for ay . The terms in the first hne of 
Eq. (jUUI) are NLO and we have to consider double insertions of such terms. All 
other terms are NNLO. Insertions of the potentials with D^^l, -D^?^, D^^'^ and 

(2) 

D2 s result in divergences. Thus, the corresponding coefficients have to be known 

in d dimensions. The insertions with C4, D^l and D^l can be related to other 
insertions, discussed below, using the (i-dimensional equation 

(^-e]g,{p,p';E)= (61) 
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For the terms in the second line of Eq. (j6(jp we get 
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where we defined q^j = p^ — pj, qij = |pj — Pj| and 

SVeom = -C4^(27r)'^5('^)(q23)-(2c4a. + 4Dfi)7rC^-^ (63) 
C^a.vr^ / C4 , ^^(2)\ e^^- T^^ - e)r(i + e) 



mqlt' V"'^ 2 + ^^^'V /i^^e ;r3/2r(l - 2e) 
We thus need the following insertions: 
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All these insertions have been computed in Ref. For the reader's convenience 
we present the results here. We write them in terms of 
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n2 = 4i^5(l,l,l,l,l;2,2,2,l-A;l) 
and the argument of the functions [iplx] = dlogr[x]/dx) and its derivatives is 



always understood to be (1 — A) = l — {Cpas)/{2y—E/m) unless stated otherwise. 
We start by writing the Green function before MS subtraction 
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where Q*^ denotes the term of order e of the Green function. This term is not 
explicitly known, but it is not needed for the calculation of the cross section, 
even though it results in finite terms in some of the insertions below, because it 
cancels in the total sum. 

The insertions of the higher-order corrections to the Coulomb potential with 
terms of the form log-'[g^/ n1]/q'^ are computed by taking derivatives with respect 
to K at K = of JUdpi G'c(pi,P2) (g23/Ats)"^"^''Gc(P3,P4)- Taking the first 
derivative yields the insertion of a single logarithm. We get 
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The insertion of a logarithm squared is obtained by taking the second derivative. 



2e3 



+ ^ + - + 



71 



12 e 



(76) 



27 



- 16A^'2 -64A(^ + i -2- ^)^'- 128^ 

- 32(^ - a) ^" + yA V'"' - 1287^2 - 64(l - i>) Hi 
+ 16 5i + 52C(3) + 196 + — -, 

where 

^1 = - E (1 + n)^r[[V^^L A] - + ^ - A] + (1 + + ^]) ■ 

(77) 

Double insertions of Coulomb potentials are computed in a similar way. The 
results read 
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Inserting 1 results in the square of the Green function 

Ii = [G{0,0;E)]\ (82) 

Note that Ii obtains a finite contribution due to The insertion due to the 5 
function can be computed directly and takes a very simple form 
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while the insertion due to the non-analytic potential is more complicated and, 
written in terms of 

£ = £ - 3 + log 2 (84) 

reads 
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Finally, the insertion due to the operator, Eq. (fTUI) . can be related to X„a and 
Xi with the help of Eq. (fTH|) and we get 

J^. ^ |j|dp^G,(pi,p2) |^(^pI-^^ -lja(P3,P4) (86) 
^2 m^a? / 1 „ 1 1 



- V + 3 



2(47r)2 V4e 2A2 A 

We are finally in a position to write down the expression for the Green function 
G™i'-'- = (87) 



2m 



— I 



2 ( Cpasir"^ 
2m 



, ,... r2(i-e)r(i + e) 

7r3/2r(l-2e) 

We mention once more that we use strictly expanded results. Thus the terms 
beyond NNLL that are present in Eq. ()87|) are dropped. 



B RGI potentials 

For ease of reference we explicitly display the RGI potentials not shown in the 
main body of the paper (we remind that ^us = f^'i/f^h)'- 
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where {z 



CoifJ's) 



A ^3 



6/3o 



«s(/^s)log 



as (/is) ll + ^§^log 



3/3o 



Qis(/^s), 



as{iJis)c\{iJis) - 



3 

~ 1 - l/{27r)aM log( J^), Po = - iTpUf) 

2z'^'' - 1 , 

9Ca + STirn/ \ AC A + 4rFn / ^ 2{Ca- 2T^mf ) 

-7C^ + 32CaCf - ACATpUf + 32Cj7Tpn/ 4rp.^^/3_2(g^/3 
A{CA + TFnf){2TFnf-CA) 



STFUf 
' 9Ca 



or, /20 32Cj.\ 



1 — Z 6 
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{2Cf - 3Ca) -tC^sM 



27CI 



TT 



tts(/i^ 



9CA + 8T^n//5o ^ 
+ 16Cf - STFUf 



5Ca + ATfUj (3o (^p^-2Ca 
ACA + ATFnf(5o-2CA 

- l) 



2{CA-2TFnf) 
-7Cl + 32CaCf - ACATFTif + 32CFTFnf 
A{CA + TFnf)i2TFnf-CA) 



X 



3/3o 



3Po + ^TFnf-2CA 
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^^Po+4TFnf/3-2CA/3 _ l 



' 9Ca 



1 



6/3o 



r- 



6/3o - 13Ca 
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